Influence of oxygen vacancy on the electronic structure of Hf02 film 
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We investigated the unoccupied part of the electronic structure of the oxygen-deficient hafnium 
oxide (HfO^i.g) using soft x-ray absorption spectroscopy at O -R" and Hf edges. Band-tail states 
beneath the unoccupied Hf 5d band are observed in the O A"-edge spectra; combined with ultraviolet 
photoemission spectrum, this indicates the non-negligible occupation of Hf 5d state. However, Hf 
A's-edge magnetic circular dichroism spectrum reveals the absence of a long-range ferromagnetic 
spin order in the oxide. Thus the small amount of d electron gained by the vacancy formation does 
not show inter-site correlation, contrary to a recent report [M. Venkatesan et al, Nature 430, 630 
(2004)]. 

PACS numbers: 78.70.Dm, 73.20.Hb, 73.61. Ng, 75.70.-i 
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Oxygen vacancy in a large-bandgap oxide changes the 
electrical characteristics by forming defect states inside 
the bandgap. Depending on the position and density of 
the defect energy levels, the chemical potential of the 
oxide - frequently identified with the Fermi level (Ep) 
- changes as in the case of doping with other atomic 
species. Hence, the characterization of the defect level 
is a prerequisite for understanding the electron transfer 
across the oxide-semiconductor interface, because the en- 
ergy barrier height at the interface depends on Ep of 
the oxideii^. With regard to hafnium oxide, the en- 
ergy position of this defect level has been estimated 
theoretically^'— and experimentallj*^ii to be above the 
middle of the bandgap, so Ep is shifted upward in the 
presence of oxygen vacancies. Since Ep is determined 
from the charge neutrality condition, it is largely influ- 
enced by the detailed density of states (DOS) near the 
bandgap^; in this case, the intensity of the features near 
the conduction band (CB) is expected to be higher than 
those near the valence band (VB). 

Recent density functional theoretical (DFT) 
studies2ii^i^i2ii2iiiS have revealed that the defect state 
is the localized Hf d state. The nonzero occupation in 
Hf d shell is very natural, because Hf atom no longer 
donates electron to the adjacent vacant sites. If an inter- 
site correlation existed between the d-electrons, even a 
long-range correlation in the electronic configuration of 
each site might emerge. Surprisingly, room-temperature 
long range ferromagnetism (FM) was recently reported 
in an oxygen-deficient Hf oxide filroiiii^ii^, although 
an experimental counter-evidence shows the absence of 
FMi^. Since the measured magnetization is the macro- 
scopic quantity, it might significantly depend on the 
purity of the sample. Furthermore, recent DFT studies 
have revealed that FM can be induced by a hafnium 
vacancy rather than an oxygen vacancy, while the oxygen 
vacancy can at the most yield paramagnetis m^^i^^i^^ . 
Thus, it is important to investigate the possible long- 
range spin order in the oxygen-deficient HfOa:<2 from 
a microscopic point of view, in that the microscopic 



quantities such as the spin or angular momentum drive 
the macroscopic properties, i.e., the magnetization. 
Since these microscopic values are determined from the 
electronic structure, the study on the local electron 
configurations of each atomic species will provide the 
information that is not influenced by a macroscopic 
deformation (domain interaction, dislocation, etc.) or 
an intervention of impurities. 

Therefore, in this paper, we focus on electron redis- 
tribution and the possible long-range spin order in the 
presence of an oxygen vacancy in Hf-oxide by studying 
its electronic structure. To investigate the unoccupied 
electronic structure of Hf d states, we performed x-ray ab- 
sorption spectroscopy (XAS) study on the HfOx<2 film. 
Since XAS reflects a transition probability from a spa- 
tially localized core-level to the unoccupied levels, it is 
sensitive to the unoccupied states originating from the 
atoms near the photon-absorbing atom. Further, owing 
to the "locality" that allows transitions irrespective of 
momentum transfer, it reflects the momentum- averaged 
electronic structure that would be equivalent to the par- 
tial density of states (PDOS) in the absence of an elec- 
tron correlation. We can choose the photon absorber by 
tuning the photon energy range to any of the core-level 
energies of each atomic species. Using the photon energy 
range ^ 530 eV, we can measure O K-edge (Is ^ 2p) 
absorption to investigate the O 2p PDOS. Since parts of 
O 2p wavefunctions are hybridized with Hf 5d wavefunc- 
tions, O if -edge XAS also shows the unoccupied Hf 5d en- 
ergy levels that are hybridized with the O wavefunctions. 
This implies that an Hf d-electron state that stays near 
the vacant site is hardly observed in O K-edge XAS. The 
situation is described in Fig. [T] The dashed circle near 
an oxygen atom roughly indicates the probing area of O 
K-edge XAS. However, at Hf TVs-edge (4^3/3, ~ 380 
eV) , we can detect all of the Hf 5d wavefunctions whether 
the d-electron state is localized near the vacancy or hy- 
bridized with the surrounding oxygen atoms (see Figure 
[T]). Although the obtained spectra become complex due 
to the on-site core-hole effect, this measurement directly 



shows the d-occupation at the Hf site. 

The O K-edge and Hf iVs-edge XASs of the oxygt 
deficient Hf oxide were performed at the 2A EPU6 beai 
hne of the Pohang Light Source (PLS). In order to c 
termine the existence of a long-range FM spin order, ■ 
also performed magnetic circular dichroism (MCD) me 
surements. The MCD spectrum is the difference betwe 
two XAS spectra measured using external magnetic fiel 
of mutually opposite directions. Here, a circularly poL 
izcd photon should be used, since the relative oricni 
tion of its chirality with respect to the direction of t 
spin moment determines the preference of pumped spii 
To align the spin direction, the magnetic field excee 
ing the coercivity field of the system should be applie 
The oxygen-deficient film of ~20 A was fabricated onto „ 
HF-etched Si(lOO) substrate with a 3 Hz pulsed Nd:YAG 
laser at a substrate temperature of 700 °C. The base 
pressure of the growth chamber was 1x10^^ Torr and no 
oxygen flow was supplied during the fllm growth. The 
oxygen concentration deduced from the intensity ratio of 
O Is/Hf 4/ core-levels in the in situ x-ray photoemission 
spectroscopy (XPS) spectra was approximately 1.8 (i.e. 
10% oxygen deficiency). An approximate estimation of 
the Hf valence leads to '-~^Hf^-^+, i.e., a composite of 60% 
of Hf4+ and 40% UP+ . 

Figure [2] shows the O A'-cdge XAS spectrum with cir- 
cularly polarized photons for HfO^i.s. The overall fea- 
tures of the spectrum are similar to those of the previ- 
ous experiments for HfO ^^^i^^'^" as well as the theoreti- 
cal band calculations for the unoccupied state9i2i2ii22i2^. 
The spectrum was analyzed using a model cluster calcu- 
lation with full atomic multiplets^l, assuming the com- 
posite [^"(60%)©^^ (40%)] in a perfect monoclinic crystal 
structure. As already noted, XAS is relevant to the in- 
teractions between the near atoms only; therefore, the 
single cluster model with an Hf atom and its nearest 
oxygen atoms is sufficient to simulate the XAS spectra. 
Though the Hf oxide film is in an amorphous phase, the 
local bonds can be assumed to be similar to those of the 
bulk monoclinic structure^^. The crystal field (CF) split 
energy ratios were calculated based on the bulk crystal 
structure, and the overall CF splitting energy was used as 
a variable scaling parameter to fit the experimental XAS 
spectra. The one-electron energy level for Hf 5dP 5d^ 
(empty bars in Fig. [5]) is analyzed to be composed of five 
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non-degenerate states: ^D{xy/x — y ), S{3z^ 



'^P[xz/yz), and two mixed states in the increasing order 
of energy. (Here the coordinate is taken as in Fig. O) 
The intensities of the five d levels are multiplied by the 
Hf 5(i— O 2p hybridization strength using the Harrison's 
rule (exponent -7 for different interatomic distances)^^ 
for the local structure as shown in the figure. The hy- 
bridization strength was found to exhibit an anisotropy; 
the in-planar state has much weaker hybridization 
strength than the other levels (~ 20% of "^S and ^P's). 
The Sd""^ — > 5(i^ simulation (filled bars) shows the similar 
results, except an overall broadening of features due to 
the atomic multiplets in the final state^^. The spectra 
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FIG. 1: Schematic diagram of the atomic distribution of 
Hf02,<2/Si. For simplicity, all the bonds are assumed to be 
four-fold. Lobes indicate the inter-site wavefunction overlaps. 
The shape and darkness in each lobe represent the electron 
density. The dashed circles indicate the spatial extent to 
which each XAS spectrum could detect the electronic struc- 
ture. In O K-edge XAS, it is impossible to observe the local- 
ized Hf d-electron state near the oxygen vacancy (denoted as 
"V"), if it were not for the delocalization of the wavefunctions, 
in contrast to the Hf A'3-edge XAS. 
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FIG. 2: Unoccupied DOS of HfO^i.8 probed by O K-eAge 
XAS with circularly polarized photons. The results of the 
atomic single cluster model calculations [d^{d^) —* d^{d?)\ are 
also shown as empty (filled) bars. The cluster considered in 
this study is shown in the inset. For comparison, the Hf 5d 
PDOS of Hf02 taken from Ref.— is shown at the bottom of 
figure. 



simulated with the overall CF energy of approximately 4 
eV, agree with the experimental spectrum and with the 
calculated Hf 5d PDOSi shown in Fig. [TJ 

A slight difference in the XAS spectra of the oxygen- 
deficient film from the stoichiometric film can be seen 
in the features near the bandgap as shown in Fig. [31 
For comparison, we appended the XAS spectrum of the 
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FIG. 3: Combined (normalized) UPS and XAS spectra of 
HfOo^i.g and stoichiometric Hf02, for DOS near the bandgap. 
The UPS spectrum for HfOi.s is taken from Refi^. Here, the 
position of the chemical potentials (zero energy in the figure) 
for both UPS and XAS spectra were arbitrarily chosen to be 
approximately in the middle of the bandgap, preserving the 
bandgap of bulk Hf02(~5.7eV). Dotted lines show the extrap- 
olations to obtain the approximate energy positions (arrows) 
of the CB and VB edges, and the shaded area indicates the 
contributions of the Si substrate or the interface state. 
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(a) Experiment: Hf A^j-edge XAS 
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FIG. 4; (a) Hf A'^s-edge XAS with circularly polarized pho- 
tons. The difference spectrum or MCD has no feature, 
indicating no long-range FM spin-order. The results of 
the atomic cluster model calculations for the composite of 
d°(60%)ed^(40%) are also shown in (b). 



stoichiometric Hf02 film measured under the same condi- 
tions. Both spectra were normalized approximately pre- 
serving the total absorption intensity. The spectrum of 
oxygen-deficient film appears broader than that of the 
stoichiometric one. We checked the broadness for many 
measurement points on each sample and found the con- 
sistency for each sample. Thus, the broadness can be 



regarded as an intrinsic property of the oxygen-deficient 
film. Though no prominent gap state is observed in the 
XAS spectrum of the oxygen-deficient film, the CB edge 
extends toward the bandgap as compared to the stoi- 
chiometric one. This can be interpreted as the partial 
delocalization of the defect state; that is, the defect state 
might be rather dispersed so that it can also be detected 
at the oxygen site. Thus, the role of the oxygen vacancy 
can be regarded as donating electron back to the Hf site 
and consequently to the O site through the Hf— O hy- 
bridization. In particular, when the vacancy is at the 
Hf02/Si interface, it is easily substituted by the Si atom 
to form an O— Hf— Si bonding. This could also slightly 
enhance the Hf d occupation, although its amount is re- 
duced due to the electron transfer to the Si site^*. How- 
ever, the CB tail cannot be attributed to the formation 
of a direct O— Si bonding. This is because the bandgap 
of SiOz (~9 eV) is much larger than of HfOa (~5.7 eV); 
therefore, no structure can exist near the bandgap. We 
can expect that the bandgaps of any types of the their 
intermediate state, i.e. the silicate (RixSiy04) will be 
between these two values^^i^S, resulting in the absence of 
states inside the Hf02 bandgap. Therefore, it is evident 
that the lack of oxygen atom in the oxide or at the in- 
terface is responsible for the occurrence of the CB tail 
stated. 

If degree of amorphousness in the Hf-oxide is somehow 
enhanced, this could also extend the CB edge due to the 
broadening of the features^^. In this case, the tailing 
should occur in both the VB and CB edges. In the fol- 
lowing paragraph, however, we will show that the tailing 
does not occur in the VB edge; this manifest that the CB 
tail is truly caused by the physics of defects, not by the 
uncertainty in the lattice dynamics^. 

To verify the existence of the band-tails in the occupied 
DOS, we also performed ultraviolet photoemission spec- 
troscopy (UPS) using He— I {hv = 21.2eV) light source 
for both the films. (The UPS spectrum for HfO^i.s is 
taken from Rcf.''^.) The UPS spectra are also normal- 
ized with respect to the total intensity. Here, the energy 
zero is arbitrarily chosen to approximately sustain the 
energy difference between the lowest XAS edge (~-|-2.5 
eV) and the highest UPS main feature (<-3 eV) to be 
the known bandgap of Hf02(~5.7 eV; Ref,^). The UPS 
spectrum of HfO^i g film is shifted by -1-0.3 eV in order 
to compensate for the Ep shift in the oxide, which was 
obtained from the cutoff kinetic energies of the secondary 
electrons in the spectra (not shown here)^; this Ep shift 
induces the higher binding energy (BE) shift in the pho- 
toemission spectra, because the BE is, by definition, the 
energy difference between the Fermi level and the core 
level, while XAS is independent of the Ep position be- 
cause the energy difference between the initial core level 
and the final conduction level is independent of the Fermi 
level. The large and broad features around the relative 
energy of -8 ~ -4 eV correspond to the main O 2p oc- 
cupied state, while the tail toward the bandgap (shaded 
area at around -3 eV in Fig. [3]) probably originates from 
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the underlying Si wafer'^, because this VB tail is absent 
in the UPS spectra of thicker (> 100 A) HfOa film (not 
shown here). A slight increase in the tail in the oxygen 
deficient film might be a signature of the metallic Hf sili- 
cide (Hf— Si bond) at the interfacc'^^ . However, the VB 
edge of the oxygen-deficient film does not extend inside 
the bandgap as compared with that of the stoichiometric 
sample. Though it is difficult to unambiguously analyze 
these band-tail states in both the occupied and unoccu- 
pied DOSs, we can observe that only the unoccupied DOS 
near the bandgap tends to extend inside the bandgap in 
the presence of oxygen vacanc y^°i^^ . This result should 
be correlated with the lifting of Fermi level because the 
virtual transition to CB becomes easier. 

Up to this point, we have shown experimental evidence 
that the oxygen vacancy lowers the CB edge and increases 
the effective Fermi level. The detailed mechanism ex- 
plaining how the vacancy formation lowers the CB edge 
could be revealed by a further consideration of the atomic 
rearrangements or the relaxation following the vacancy 
formation. For example, the disorder in the positions of 
the vacant sites can play a role in broadening the CB DOS 
and consequently in reducing the bandgap. Further, the 
energy position of the defect level is found to depend on 
the atomic arrangements after the vacancy formation^. 
The reduction in the mean Hf— O bond length (^-0.1 
A) and in the number of bonds, in the presence of oxy- 
gen vacancy^^, will influence the detailed balance in the 
charge dynamics between the two atomic species. The 
CB lowering can also be correlated with the broadening 
of features influenced by the valence conflgurations (the 
number of electrons) of the remaining atoms, as shown 
in the simulation results in Fig. O Further investigation 
is required to completely understand the mechanism. 

Now we check if there be a long-range FM spin order. 
The Hf TVs-edge XAS spectra obtained with circularly 
polarized photons are shown in Fig. SJa). The external 
magnetic field was maintained perpendicular to the sam- 
ple plane with H — ±2500 Oe, which is much larger than 
the reported ferromagnetic coercivity field of the Hf02 
film^^. The difference spectrum, i.e., the MCD spectrum 
exhibits no clear feature and remains zero within the ex- 
perimental uncertainty. This indicates that there is no Hf 
d-electron-related FM spin order. One would argue that 
the density of the oxygen vacancy is too small to detect 
FM; however, the value of 10% is larger than the dopant 
density in conventional diluted magnetic semiconductors 
(DMS)^. 

To analyze the spectral features, we appended in Fig. 



Illjb) the simulation result obtained by the full atomic 
multiplet calculation for Hf iVs-edge XAS 5(i) 
with the parameters similar to that of O i^-edge XAS, 
except the 4p — 5d Coulomb interactions'^^. As was in the 
O K-edge XAS simulation, the simulated spectra are ob- 
tained by the sum of two transitions: 5(f — > Ap^5d^ and 
4p'^5(P. The features of the latter are shifted by 
U — — —2 eV''^ with respect to those of the former. 
The result of the calculation is in an excellent agreement 
with the experimental spectra; therefore, it provides def- 
inite peak assignments; the Hf 4p3/2~*5(i transition is 
the largest feature, while the other features in the higher 
photon energies are attributed to the Hf 4p3/2— >6sp tran- 
sition. Contrary to 2p — > 3d transition (as in the case of 
transition metal ions), the feature of Hf Ap 5d tran- 
sition is not much stronger than ip —f 6sp, because of 
its lower photo-absorption cross-sections. The features 
of Hf 6sp for both the edges are denoted by T in Figs. 
[2]and|31 The peak positions almost coincide, confirming 
our peak assignments. The simulated MCD signal for Hf 
5d [Fig. SKb)] is small due to the small anisotropy in the 
unoccupied electron states {d^~^). Although the sensi- 
tivity of the MCD measurement on magnetism is low, no 
correlation was observed between the experimental and 
simulated MCD spectra. This confirms that the small 
modulation in the experimental spectrum is caused by 
the noise, not the possible long-range magnetism. 

Therefore, we can conclude that there is no preference 
in spin selection in the electron occupation to either Hf 
5d or 6sp; this evidently proves "no ferromagnetic spin- 
order" in the case of oxygen vacancy. However, the elec- 
tron occupation number at an Hf site should be nonzero; 
otherwise, the electronic configurations of the next near- 
est O sites, which was shown in O i^-edge XAS spectra 
(Fig. [3]), should never be influenced by the presence of 
oxygen vacancy. Thus, the absence of FM spin-order can 
be explained by the absence of long-range correlation be- 
tween the d-electrons in the neighboring Hf sites, rather 
than the absence of d-electron itself. 

We conclude our arguments with two remarks: the 
oxygen vacancy in Hf oxide slightly enhances the CB 
DOS. However, it docs not involve the long-range ferro- 
magnetic spin order in the oxide because of the absence 
of inter-site electron correlation. 
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